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The coupling between flame surface 
dynamics and species mass conservation 
in premixed turbulent combustion 

By A. Trouv^ 1 , D. Veynante 2 , K. N. C. Bray 3 AND T. Mantel 4 


Current flamelet models based on a description of the flame surface dynamics 
require the closure of two inter-related equations: a transport equation for the mean 
reaction progress variable, c, and a transport equation for the flame surface density, 
E. The coupling between these two equations is investigated using direct numerical 
simulations (DNS) with emphasis on the correlation between the turbulent fluxes 
of c, pu n c n , and E, (u") 5 E. Two different DNS databases are used in the present 
work: a database developed at CTR by A. Trouve and a database developed by 
C. J. Rutland using a different code. Both databases correspond to statistically 
one-dimensional premixed flames in isotropic turbulent flow. The run parameters, 
however, are significantly different, and the two databases correspond to different 
combustion regimes. It is found that in all simulated flames, the correlation between 
pu"c" and ( u n ) S E is always strong. The sign, however, of the turbulent flux of c or E 
with respect to the mean gradients, ffc/dx or dYi/dx, is case- dependent. The CTR 
database is found to exhibit gradient turbulent transport of c and E, whereas the 
Rutland DNS features counter-gradient diffusion. The two databases are analyzed 
and compared using various tools (a local analysis of the flow field near the flame, a 
classical analysis of the conservation equation for u n c n , and a thin flame theoretical 
analysis). A mechanism is then proposed to explain the discrepancies between the 
two databases and a preliminary simple criterion is derived to predict the occurrence 
of gradient /counter-gradient turbulent diffusion. 


1. Flame surface density and species mass evolution equations 

The objective of theoretical descriptions of turbulent reacting flows is to provide 
tractable expressions for unclosed terms appearing in the conservation equations 
for mass, momentum, and energy. In the following, we focus our attention on 
mass conservation. In the classical theory of turbulent premixed flames, under the 
assumption of simple chemistry, the mass fractions of the reactive species are all 
linearly related and may be expressed in terms of a single reduced mass fraction 
called the reaction progress variable, c= 1 - ( Yr/Yr >u ), where Yr is the fuel mass 
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fraction and Yr iU its value in the unburnt gas. The ensemble-averaged species mass 
balance is then written as: 

%jjp- + V.(pUc) + V.^P 7 ? 7 ) = V.(^PVP) + u r /Yr iU , (1) 

at 

where the tilde superscript denotes a Favre (density- weighted) average Q = pQ /p, 
the double prime symbol denotes the instantaneous deviation from the Favre aver- 
age, Q” = Q — Q, D is the fuel mass diffusivity, and Cjr represents the mass of fuel 
consumed by chemical reaction, per unit time and per unit volume. 

In the flamelet theory for turbulent premixed combustion, the reaction zone is 
assumed to be a thin surface separating fresh and burnt gases. The local reaction 
rate may be expressed in terms of the local flame surface-to- volume ratio, E', and 
the ensemble- averaged reaction rate may be expressed in terms of tjie mean flame 
surface- to- volume ratio, also called the flame surface density, E = E* (Bray 1980, 
Williams 1985, Peters 1986): 


wfi = p%YR t% (Sc)s £, ( 2 ) 

where p u is the density in the unburnt gas and (Sc) $ ls fhe mean fuel consumption 
speed. The mean consumption speed accounts for local variations of the reaction 
rate along the flame surface while the flame surface density, E, characterizes the 
flame wrinkling due to the turbulent motions. For flames with Lewis numbers close 
to unity, {Sc)s remains close to the laminar burning velocity si and to first order 
the mean reaction rate is proportional to E (Haworth & Poinsot 1992; Rutland & 
Trouve 1993; Trouve & Poinsot 1994). The c-equation is then re-written as: 

- + V.(pUc) + V.(^ I 77 ? 7 ) - V.(^0VP) + P uS L X. (3) 

at 

In (3), the contribution of molecular diffusion is usually neglected for high Re ynolds 
number flows and closure is only required to describe the turbulent flux of c, pu"c' / 
and the flame surface density, E. A variety of modeling choices may be made f or 
those terms ranging from standard gradient transport approximatio ns for pu ft c n or 
simple algebraic closures for E, to full transport equations for both pu ,, c n and E. 

For instance, in the Coherent Flame Model (Marble & Broadwell 1977; Darabiha 
et al 1987; Maistret et al 1989; Candel et al 1990) the flame surface density 
is obtained via a modeled formulation of an exact evolution equation called the 
E-equation (Pope 1988; Candel & Poinsot 1990): 

^ + V.UE + V.(u") s E + V.(um) s E = (k)s E, (4) 

at 

where w is the flame front propagation speed, n is the flame normal vector pointing 
towards the fresh gases, k is the flame stretch, and ( )s denotes a fla me su rface 
mean defined as an area- weighted ensemble- average (Pope 1988), (Q)s = Q^ f / S' = 
QE'/E. The three convective terms on the left-hand side of (4) are transport terms 
that correspond respectively to convection by the mean flow, turbulent diffusion, 
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and flame propagation. The term on the right-hand side of the equation is the 
source/sink term for flame surface density and accounts for production of flame 
surface area due to hydrodynamic straining and dissipation due to combined effects 
of flame propagation and flame surface curvature (Trouve & Poinsot 1994). Different 
closure assumptions are required in the E-equation, in particular to calculate: the 
turbulent diffusion velocity, (u") s , the transport due to flame propagation, (u;n) 5 , 
and the turbulent flame stretch, (k)s (see Duclos et al. (1993) for a critical review 
of the different formulations of the modeled S-equation that can be found in the 
literature). 

While the need to model E through a transport equation has become increas- 
ingly recognized in recent years, the modeling of the turbulent flux of c, pu "c", 
remains in comparison somewhat controversial. Using a standard gradient diffusion 
approximation (GD), one may write: 


p “> c " = <5) 
where p t is the turbulent viscosity. 

Premixed flames, however, are known to exhibit counter-gradient diffusion (CGD) 
of mass reactant (Libby & Bray 1981; Bray et al. 1981). CGD is related to the 
differential effect of mean pressure gradients on cold, heavy reactants and hot, light 
products. This effect may be shown very simply under the classical Bray-Moss- 
Libby assumption of fresh reactants (c = 0) and fully burnt products (c = 1) 
separated by a thin flame sheet, leading to a bimodal probability density function 
of c. The x-coinponent of the turbulent flux may then be written as (Bray 1980): 

= pJ"?' = pc( 1 - c) (U b - U u ) , (6) 

where U b and U u are respectively the x-component of the conditional mean veloc- 
ity within burnt and unburnt gases. Thermal expansion and the associated flow 
acceleration through the flame will tend to make Ub greater than I7 U , thereby pro- 
moting counter-gradient turbulent diffusion of c {pu u c n /Sfcjdx > 0, contrary to the 
predictions from Eq. (5)). Note that counter- gradient diffusion has been observed 
in a number of experiments (Moss 1980; Shepherd et al 1982; Cheng & Shepherd 
1991; Armstrong &: Bray 1992). 

While CGD has been extensively studied in the past ten years, the implications 
of this non-gradient transport phenomenon on the distribution of flame surface 
densities and the overall mean reaction rate remains unknown. In fact, a variety 
of conflicting modeling assumptions are still used in the current literature ranging 
from assumptions that simply neglect CGD and use standard gradient transport 
approximations as displayed in (5) (Darabiha et al 1987; Maistret et al 1989; 
Candel et al 1990; Cant et al 1990) to formulations where such approximations 
are carefully avoided and where closure is achieved by writing a transport equation 
for p\i u c n (for instance in the Bray-Moss- Libby model: Bray 1980; Bray et al 1989; 
Bray 1990). 



98 


A. Trouve, D . Vcynante , K. N. C. Bray & T. Mantel 


It is important to realize that the closure assumptions required in the E-equation 
cannot be made independently of those required in the equation for c. This may be 
best understood by considering the interdependence of E and c. As seen in Eq. (3), 
E is a source term for c and the mean reaction progress variable will depend strongly 
on the distribution of flame surface densities. In addition, following Pope (1988), 
E may be expressed in terms of statistical properties of the c-fiekh For instance, 
the flame surface density is expressed as the product of the expected value for 
the magnitude of the gradient of c, conditioned on the flame surface, times the 
probability of being on that surface: 

£ = <|Vc| | c = c f ) p(c f ), (7) 

where the flame is viewed as a surface-contour, c = cj, and where p{cf) is the 
probability of c = c/ . 

This theoretical coupling between E and c must have implications for models. 
Some of those implications are already well-known. For instance, in current models 
of the E-equation, while the turbulent transport, {u") s , is described exclusively in 
terms of statistics of the turbulent flow field, the realizability of the combustion 
model is ensured by using a modeled expression for the flame stretch, (k)s , that is 
c-dependent. For instance, the closure model for (k)s leads to a vanishing E when 
the reaction reaches completion and c approaches unity. 

The objective of the present work is to study further the exact implications of 
the coupling between the equations for E and c. More specifically, direct numerical 
simulation (DNS) is used to determine: (1) if and under what conditions gradient- 
diffusion (GD) or counter-gradient diffusion (CGD) is observed in the c-equation; 
(2) if there is any correlation between the turbulent flux of c and the turbulent 
flux of E, and in particular whether CGD (GD) in the c-equation occurs with CGD 
(GD) in the E-equation. 

2. Two different DNS databases 

In the present work, two different direct numerical simulation databases are used. 
The first one, referenced as the CTR database, was developed at CTR by A. Trouve 
and was previously used to study the statistics of the turbulent flame stretch, (k) s , 
across the turbulent flame brush (Trouve & Poinsot 1994). The CTR database was 
developed using a fully compressible, 3D, Navier-Stokes solver with one step irre- 
versible chemistry. The numerical configuration corresponds to a premixed flame 
propagating in decaying isotropic turbulence. The ratio between the turbulent rms 
velocity and the laminar flame speed, u'/si, is set to an initial value of 10 and 
decreases to about 3 after four turbulent eddy turnover times. The second database 
was recently developed by C. J. Rutland (see Rutland fc Cant 1994) using a low 
Mach number, 3D code. The Rutland code assumes a constant viscosity, whereas 
the CTR code features temperature- dependent transport coefficients. The numer- 
ical configuration corresponds to a premixed flame propagating in approximately 
stationary, weak turbulence, u ' /sl ~ 1. In the Rutland configuration, like in the 
CTR one, the simulated turbulent flame is statistically ID. Comparisons between 
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FIGURE 1. Comparison of the CTR and the Rutland DNS. The mean progress 
variable c is plotted as a function of x-location along the direction of mean prop- 
agation. The comparison is performed at a time selected so that the turbulent 
flame speed is approximately the same in both cases. Length scales are made non- 

dimensional by the laminar flame thickness, Si . : CTR database; : 

Rutland database. 



x location 


FIGURE 2. Comparison of the CTR and the Rutland DNS. The flame surface den- 
sity S is plotted as a function of x-location along the direction of mean propagation. 
: CTR database; : Rutland database. 
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the two databases are performed assuming that the different computations corre- 
spond to the same laminar flame (same laminar flame thickness, 61 , same laminar 
flame speed, sl, and same molecular transport coefficients). 

While the two databases feature turbulent flames with similar values of the tur- 
bulent flame speed, St, they also correspond to turbulent flames with significantly 
different structure. Fig. 1 shows two instantaneous profiles of the Favre-averaged 
progress variable c as a function of r-location along the direction of mean propa- 
gation. Fig. 2 presents a similar comparison for the flame surface density profiles. 
It is seen that the flame brush is about three times thicker in Rutland’s database 
compared to the CTR case. Also, the flame front wrinkling, as measured by the 
magnitude of the flame surface density, is much greater in the CTR database. As 
seen in Fig. 3, the two DNS lie at different locations in the classical turbulent 
combustion diagram due to Borghi and Barrere (see Borghi 1985) and therefore 
correspond to different turbulent combustion regimes. The CTR database corre- 
sponds to flames that are more turbulent and feature smaller length scales. The box 
plotted in Fig. 3 corresponds to a zone where counter-gradient diffusion has been 
experimentally observed. Rutland’s simulation lies in this zone, the CTR database 
lies outside. Note also that according to the Klimov- Williams criterion, the CTR 
database corresponds to non-flamelet combustion. Recent work by Poinsot e< al. 
(1991) has shown that the domain of flamelet combustion is in fact significantly 
larger and the CTR database is within that domain according to the Poinsot et al. 
criterion (see Fig. 3). 

The most interesting result is that the two databases display striking differences 
in their turbulent transport properties. Fig. 4 shows the turbulent flux pu"c" as 
a function of the mean progress variable c. Rutland’s database exhibits counter- 
gradient turbulent diffusion of c (pu"c" > 0), whereas gradient diffusion transport 
is found in the CTR database ( pu"c " < 0). These differences can be related to dif- 
ferences between the two codes (compressible/incompressible and variable/constant 
molecular transport coefficients), differences due to different initial and boundary 
conditions, or differences due to different values of the run parameters corresponding 
to different combustion regimes. 

Figs. 5a and 5b present the spatial variations of different relevant mean flow 
velocities, U, U u , V b, and (u) s , across the turbulent flame brush. In the Rut- 
land simulation, the mean velocity within the products, U b, is always greater than 
the mean velocity within the reactants, U u , which according to (6) corresponds to 
counter-gradient turbulent diffusion of c. In the CTR database, however, U u is 
greater than Ub, which again according to (6) corresponds to gradient turbulent 
diffusion. This last result may seem surprising since it is expected that the thermal 
expansion will induce a burnt gas velocity larger than the fresh gas velocity. It is 
worth emphasizing that conditional velocities may be difficult to interpret because 
the sampling is quite different for U u and for Ub- For example, at the leading edge 
of the turbulent flame brush, the flow field corresponds mainly _to values of the 
progress variable c = 0, and the fresh gas conditional velocity, U u , is computed 
from a large number of samples. At that location, a few pockets have a value of 
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FIGURE 3. Premixed turbulent combustion diagram showing the different com- 
bustion regimes (from Borghi 1985; and Poinsot et al 1991). The coordinates are 
the ratio of the turbulent length scale, /*, and the laminar flame thickness, <5/, and 
the ratio of the turbulent fluctuation, id, divided by the laminar flame speed, sl> 
The parameters of the two DNS databases are plotted in the diagram. 



Mean progress variable 


F IGU RE 4. Comparison of the CTR and the Rutland DNS. The turbulent c-flux, 
/5u"c", is plotted as a function of the mean reaction progress variable, c. Velocities 

are made non-dimensional by the laminar burning velocity, si. : CTR 

database; ; Rutland database 
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FIGURE 5. Mean flow velocities across the turbulent^ flame. The Favre- averaged 
velocity, U y the conditional mean velocities, U u and Ub, and the surface-averaged 
velocity, (u) s , are plotted in c-space. Velocities are made non-dimensional by the 
laminar burning velocity, s/,. Top: Rutland database. Bottom: CTR database. 

" '■■■ : Favre average; : products; : reactants; : < u >s* 
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the progress variable c = 1, and is determined from a small number of sam- 
ples. The discrepancy between Figs. 5a and 5b is one of the keys to understanding 
when gradient or counter-gradient diffusion occurs and will be examined in further 
sections. 

Different Lewis numbers Le (ratio of thermal to mass molecular diffusivities) 
have been investigatedjn the CTR database. Fig. 6 shows the spatial variations of 
the turbulent flux, u n c lf , across the turbulent flame brush at various times in the 
simulation, for two different Lewis numbers Le — 0.8 and Le = 1.0. In this figure, 
u H c n is made non-dimensional using the instantaneous rms turbulent velocity taken 
in the fresh gases. This non-dimensionalization allows separation of the decrease 
of the turbulent flux due to the decay of the turbulence from the variations due 
to^other phenomena. Two different behaviors axe observed. For Le — 1.0, the 
u"c"-profile is changing rapidly during an initial phase and remains approximately 
constant after three eddy turnover times. In contrast, for Le = 0.8 the u'V'-profile 
keeps changing in time and exhibits a continuous variation f rom negative towards 
positive values. These results Eire also shown in Fig. 7 where 'pu ,t c N is first integrated 
across the turbulent flame brush and then plotted versus time. In Fig. 7, the Le — 1 
flame reaches an asymptotic value, indicating thereby that no transition towards 
counter-gradient diffusion is to be expected and that a non-dimensional time of 4 
is sufficient for data processing of the turbulent flame brush. Gradient turbulent 
transport observed in this case does not seem to be related to effects of initial 
conditions. For Le = 0.8, there is a clear trend towards transition to counter- 
gradient turbulent transport that might be related to thermal diffusive instabilities. 

We now examine the idea that the differences in turbulent transport properties 
between the two databases are due to basic differences in the flame-flow dynamics (as 
opposed to numerical artifacts). These differences are further characterized using 
a local analysis of the flow field near the flame (section 3). Section 4 presents a 
comparison between the turbulent fluxes of c and S. A classical analysis of gradient 
versus counter-gradient turbulent diffusion based on the conservation equation for 
u”c n is presented in section 5. Finally, a theoretical analysis using a thin flame 
model is developed in section 6. 

3. Local flow structure near the flame surface 

As discussed in §2, the Le = 1 case from the CTR database and the Rutland 
simulation correspond to turbulent flames characterized by the same laminar flame 
thickness, <5^, the same laminar burning velocity, s*,, the same molecular transport 
coefficients, but embedded in different turbulent flow fields. These simulated flames 
feature different global (spatially- averaged) properties, as illustrated in Figs. 1 and 
2. A different perspective is adopted in this section where the flow velocity and 
the c-field are spatially resolved and analyzed in the vicinity of the reactive layers 
in a frame of reference attached to the flame. This frame of reference is used to 
determine, in particular, whether local flow variations occur in directions that are 
normal or tangential to the flame surface, i.e. whether local flow velocity gradients 
are aligned with local concentration gradients of reactive species. 
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FIGURE 6. Evolution of u"c" across the turbulent flame brush at different times 
in the simulations a) Le = 0.8; b) Le = 1.0. Time is made non-dimensional by the 

initial turbulent eddy turnover time. : t = 1.4; : t = 2.7; : t = 

3.6; : t = 4.5. 
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FIGURE 7 . Evolution of the x-integral of pu n c n , integrated across the turbulent 

flame brush, versus time for different Lewis numbers : Le = 0.8; : 

Le = 1.0; : Le = 1.2. 

The flame-based analysis reveals that the local reaction rate profiles across the 
flame remain approximately uniform along the flame surface and are similar in the 
CTR and Rutland simulations. In other words, the local chemical structure of 
the flame remains laminar-like, a result that is consistent with a flamelet analysis 
and shows that the chemistry of such flames with unity Lewis number is relatively 
insensitive to flow perturbations. The local flow velocity profiles, however, exhibit 
striking differences, as discussed below. 

Figs. 8a and 8b show the local flow dilatation across the flame at various locations 
along the flame surface, respectively for the Rutland and the CTR simulations. 
The dilatation of the flow is produced by both heat transfer in the flame preheat 
zone and heat release in the reaction zone. These local dilatation profiles remain 
approximately uniform along the flame surface and laminar-like in the Rutland 
simulation, whereas they exhibit more variations and significant deviations from 
the laminar case in the CTR simulation. 

In the Rutland simulation, the flow field is found to be essentially one-dimensional 
and quasi-steady close to the flame. Most of the flame acceleration occurs along the 
flame normal direction: 


V.u w V(u.n).n. (8) 

In that situation, the flow field within the flamelets is well described using the 
classical expressions that describe plane stretch-free laminar flames, where the flow 
velocity varies linearly with c in the flame normal direction: 
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a) 



FIGURE 8. A selection of local flow dilatation profiles normal to the turbulent 
flame for a) the Rutland database; b) the CTR database. In both figures, a local 
profile of the normal component of the flow velocity is also shown for comparison 
(a test of Eq. (8)). Quantities are made non-dimensional by the laminar burning 
velocity, sl, and the laminar flame thickness, 6l- 


f 
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u.n(c) ~ u.n(c') 4- r(c' — c)st, 

( 9 ) 

u.t(c) « u.t(c'), 

(10) 


where r is the heat release factor, r = T& — T tt /T u , T tt , and T b being respectively 
the temperature within fresh and burnt gases; t is a unit vector in the flame tangent 
plane. The relations (9) and (10) may be re* written as: 

(u.n(c = 0.8) - u.n(c))/rsj, + 0.8 « c, (11) 

(u.t(c) — u.t(c = 0 .8))/tsl « 0. (12) 

Expressions (11) and (12) are found to provide good descriptions of the flow 
variations within the flamelets, as seen in Fig. 9a. 

In the CTR simulation, however, the flow field is not one- dimensional and cannot 
be deduced directly from the dilatational field. Fig. 8b shows that (8) does not 
hold and Fig. 9b shows that the normal component of the flow velocity within the 
flamelets does not vary linearly with c; its gradient is not aligned with the gradient 
of c and exhibits large variations from one flame location to another. While in the 
Rutland simulation the flow field is determined by the dilatation occurring within 
the flame, the flow field in the CTR simulation appears dominated by the turbulent 
motions. These results will be used in section 6. 


4. The relation between the turbulent fluxes of c and E 

As shown in Fig. 4, the sign of the turbulent flux of c is different in the Rutland 
and CTR simulations. We now examine the turbulent flux of E, (u") 5 E. Fig. 10 
shows that the CTR database features gradient diffusion transport for the flame 
surface density, whereas the Rutland database corresponds to counter-gradient dif- 
fusion. Note that since the turbulence intensities are higher in the CTR simulations, 
the magnitude of the turbulent fluxes are also found to be higher. 

Following Bidaux and Bray (1994), an estimate of (u) s is given by: 

(u) s = U u +c*(U b ~U u ), (13) 

where c* is the progress variable level used to trace the flame front. This expression 
assumes a linear evolution of the mean velocity in the flame zone and is supported 
by Figs. 5a and 5b. Using the classical Bray-Moss-Libby (BML) relation U = 
(1 — c) U u + cU b , we can also write: 


(u")s = (*)s-U = (c*-c)(U b -U u ). 

Combining (6) and (14), we obtain: 


(u") s E 


(c*-3 

c(l-c) 


i^'c^E, 


(14) 


(15) 
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FIGURE 9. A selection of local profiles of the normal component of the flow velocity 
normal to the turbulent flame for a) the Rutland database; b) the CTR database 
(a test of Eq. (11)). The dotted line is the curve obtained for a plane, stretch-free 
laminar flame, In Fig. 9a, a selection of local profiles of the tangential component 
of the flow velocity is also shown for comparison (a test of Eq. (12)). Quantities are 
made non-dimensional by the laminar burning velocity, si y and the laminar flame 
thickness, 8i, 
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Mean progress variable 


Figure 10. Comparison of the CTR and the Rutland DNS. The turbulent E-flux, 
( u N ) 5 E, is plotted as a function of the mean reaction progress variable, c. Dashed 
lines corresponds to the estimate from Eq. (15) assuming c* = 0.5. Velocities are 
made non-dimensional by the laminar burning velocity, si. 1 : CTR database; 

: Rutland database. 

which provides a simple relation between the turbulent fluxes for E and c. This 
relation shows that (u n ) s E and u"c" have the same sign at the leading edge of 
the turbulent flame, near c = 0, and opposite signs on the burnt gas side, near 
c = 1. Gradient (respectively counter-gradient) turbulent diffusion of the mean 
progress variable, c, implies gradient (counter-gradient) turbulent diffusion of the 
flame surface density, E. The estimates obtained from Eq. (15) axe also displayed 
in Fig. 10 and are seen to provide a very good estimate of {u") s E. In any case, 
the main result is that the turbulent diffusion term in the evolution equation for E 
is strongly correlated to the turbulent diffusion term in the equation for c. In the 
following, we focus attention on the turbulent c-flux. 

5. Study of the transport equation for u"c n 

5 A Preliminaries 

To identify the physical phenomena responsible for turbulent diffusion of c, a 
classical analysis consists of examining in detail the terms of the transport equation 
for the second order moment u''c". This balance equation can be obtained in 
a classical way from the continuity, the Navier-Stokes and the progress variable 
transport equations. After some algebra, the complete equation for u'/c" may be 
expressed (Favre et al 1976; Launder 1976): 
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(I) 


(II) 


d 

dxj 

( pu'Ju'i'c ") 

— pu”u 

„ ft 

’d^~ PU] ^ 


(III) 

(IV) 

(V) 

— &P 

„dp' 

„9Jk 

, ,< dr ' k 

i — U * 

C dx, 

C dx, 

’ dx k 

+ c" 

dxk 

+ pu { w 

(VI) 

(VII) 

(VIII) 

(IX) 

(X) 


(16) 


where Sk represents the molecular diffusion flux of c and r,* is the viscous stress 
tensor. 

All the terms in (16) can be extracted from the simulations (including the un- 
steady term). For instance, Fig. 11 shows typical spatial variations of the different 
terms of (16), written for the x-component of the turbulent flux u"c", as obtained 
from the CTR database. 

Fig. 11 shows that some of the terms in (16) are dominant and will determine the 
sign of u n c n . The dominant contributions are due to the terms I, IV, VI, and VII. In 
the following section, our attention will be focused on the role played by these terms 
and their temporal evolution. Moreover, we have to notice that the imbalance term 
due to the inherent numerical approximations is one order of magnitude smaller 
than all the terms appearing in Eq. (16). This result is satisfactory and shows that 
the quality of the CTR database and that of the post-processing are sufficient to 
analyze the variations of second order moments. 


5.2 Effect of the turbulence intensity and of the mean pressure gradient 

In order to better understand the time evolution of u n c n (see Figs. 6 and 7), some 
terms in (16) are studied in more detail for two different Lewis numbers, Le = 0.8 
and Le = 1.0. These terms are: 

- the unsteady term I 

- the production due to the interaction of the mean progress variable gradient and 
the turbulent flow field (term IV in (16)). Since this term is by definition always 
positive, it is responsible for gradient type diffusion. 

- the production due to the mean pressure gradient (term VI in (16)) is supposed 
to be responsible for counter-gradient diffusion as pointed out by Libby & Bray 
(1981). These authors argue that the pressure gradient across the flame front 
preferentially accelerates the low density gases, creating a relative motion between 
fully burnt and fresh gases. 

- the correlation between the fluctuations of the progress variable and the fluctua- 
tions of the pressure gradient (term VII). Little is known about this term which 
is usually considered as a dissipation term in non-reacting flows (Launder 1976). 


Fig. 12 displays the time evolution of the terms I, IV, VI, and VII integrated 
the turbulent flame brush ( J -~pu”c n dc, J pu tt2 —dc, J 


across 
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mean progress variable 


FIGURE 11, Spatial variations of the different terms in Eq. (16) versus c, for the 
CTR case Le — 1.0 and t/ro = 4.5. The terms in (16) are made non-dimensional 
by Pu u'”/l tu . The curve with open circles correspond to the numerical errors that 
were found when closing the budget (16) and provides an estimate of the magnitude 
of the error for each term. 

For both Le = 0.8 and Le = 1.0, the unsteady term I and the dissipation term VII 
behave in the same way. At early times, the unsteady term is negative producing 
counter-gradient diffusion and becomes positive later in the simulation playing the 
opposite role. The term VII is always negative in the simulations and acts to 
promote counter-gradient diffusion. 

We now analyze the terms IV and VI which are supposed to be the key terms 
responsible for the turbulent gradient transport observed in the simulations. At 
initial times, both of these terms are positive and tend to produce gradient diffusion. 
Due to the high turbulence level (at t/ro = 1.4, u f jsi = 9.1) and the relatively 
small thickness of the flame brush, the term IV appears to be the most dominant 
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a) 


b) 


FIGURE 12. Time evolution of the terms I, IV, VI, and VII of Eq. (16) integrated 
across the turbulent flame brush, for (a) Le = 0.8 and (b) Le — 1.0. Term s I, 
IV, VI, and VII are made non-dimensional by the instantaneous value of puu'^ 2 /lt u . 
Time is made non-dimensional by the initial turbulent eddy turnover time. 

term in (16). This result is consistent with the numerical study of Masuya (1986) 
who observes the flux becoming gradient diffusion type as the turbulence intensity 
increases. 

As time proceeds, these two terms decrease continuously. While the term IV 
reaches an asymptotic value, the term VI changes sign and acts to promote counter- 
gradient diffusion. This feature is especially noticeable for the case Le = 0.8 where 
the term VI becomes the largest term in (16) after 4 turbulent eddy turnover time. 
This can be explained by the fact that for Le = 0.8 the turbulent flame speed is 
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much higher than in the case Le = 1 (TVouve & Poinsot 1994) and consequently 
provokes a higher pressure jump across the turbulent flame which is approximately 
given by A p ~ Tp u S\. 

5.8 Conclusions 

At this point, we can tentatively conclude that the different behaviors for u"c" 
observed in the CTR database for Le — 0.8 and Le = 1.0 are essentially due to a 
higher mean pressure gradient in the case Le = 0.8, which is mainly due to a higher 
turbulent flame speed. 

In a more general sense, it is well known experimentally that the turbulent flame 
speed, Sr, increases almost linearly with the turbulence intensity. Thus, the terms 
IV and VI are of the same order of magnitude and may even compensate themselves 
if the heat release parameter is sufficiently low (we can roughly estimate the ratio of 
the terms IV and VI to be of order r” 1 ). This simplistic argument is in agreement 
with the work of Masuya (1986), who shows a stronger effect of r compared to the 
effect of the ratio u' j si on the sign of u w c w . 

6. An analysis based on a thin flame model 

6.1 An expression for turbulent diffusion transport 

An analysis of the turbulent diffusion of the mean reaction progress variable may 
be conducted under the assumption of an infinitely thin flame front. Following Cant 
et al. (1990) or Trouve and Poinsot (1994), the flame is viewed as a thin surface 
propagating towards the fresh gases with a velocity X given by the sum of the fluid 
velocity and the flame propagation speed in the normal direction: X = u + u>n. It 
may be shown that: 


(JC • 

-qI = {X.n) s E = (u.n) s E + (u;) s E, (17) 

where c is the Reynolds- averaged progress variable. This equation may be re-written 
as: 


— +U.Vc=(u.n) s E + (u-) s E + U.Vc. (18) 

Using the geometrical relation Vc = — (n) 5 S, (17) becomes: 
dc ~ 

^ + U .Vc = (u".n) s E + (w) s S. (19) 

In the Bray-Moss- Libby theory, Reynolds and Favre averages of the progress 
variable c are related through the following expression: c = where pi> is 

the mass density in the burnt gas. (18) and (19) can then be recast as transport 
equations for c: 
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or 


^ + U.Vpc = p b ( u".n) s S + pb{w) s Z- 
An exact equation for the Favre- averaged progress variable is: 


( 21 ) 


^ + V.pUc = — V.pu"c" + (pui) s E, 
at 


( 22 ) 


where molecular species diffusion and chemical reaction are included in (pw) s . 

Combining (22) with (20) or (21), two expressions for the turbulent transport of 
c as a function of flame surface-averaged quantities are obtained: 


— V.pu"c" = pb( u.n) s £ + (pb(w) s ~ (pwjs ) s + V.pUc, (23) 


— V.pu"c" = pb{u" .n) S Y, + (pb(w) s - (pu») s )E + pcV.U. (24) 

These two expressions exhibit three contributions to the turbulent diffusion of c. 
The first term on the right-hand side represents the correlation between the flame 
front movement and the velocity field. The second term involves the reaction rate, 
and the third term is clearly related to the thermal expansion. One may also note 
that with a model for the reaction rate, the only remaining unclosed term is (u'\n) s . 



Mean progress variable 


FIGURE 13. Comparison of the exact turbulent diffusion term — V.pu"c" with 
its estimate from Eq. (23) in the CTR database. ■ ■" : CTR database; : 

estimation from (23). 
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Fig. 13 presents a comparison of the right-hand side and left-hand side of (23) 
in the CTR database ( Le = 1). The agreement is quite good even though the 
probability density function of c is not fully bimodal in the database as assumed to 
derive expression (23). Fig. 14 shows the correlation (u.n) 5 as a function of c for 
the two databases. This term is quite different in the two data sets. In the CTR 
database featuring a gradient turbulent transport, (u.n) 5 decreases almost linearly 
with c from positive values near c = 0 to negative ones near c = 1. The positive 
values of (u.n} 5 near c = 0 correspond to local flow velocity vectors pointing into 
the reactants. At the leading edge, the flame front is convected towards fresh gases 
by the turbulent motions. This phenomenon accounts for the negative value of Ub 
seen in Fig. 5b. The product (u.n) 5 S in (23) corresponds to a strong gradient 
trend (positive at the leading edge and negative at the trailing edge). On the other 
hand, in the Rutland database featuring a counter-gradient turbulent transport, 
(u.n) s is negative with an almost constant value. In (23), the product (u.n) 5 S is 
counter-gradient at the leading edge and becomes gradient at the trailing edge. Its 
effect is opposite to the one due to thermal expansion, \?pUc. 

Fig. 14 also shows the uncorrelated ((u) s .(n) s ) and the correlated ((u.n) s — 
(u) 5 .(n) s ) parts of the term (u.n) s . It is seen that u and n are mainly uncorrelated 
in the CTR database, whereas the correlation is not negligible in the Rutland case. 
This result is also in agreement with the local analysis developed in section 3. 

6.2 Analysis of (u.n} 5 and (u".n) 5 terms 

From the previous analysis, the terms (u.n) s and (u ;, .n) 5 are found to be im- 
portant ingredients of the turbulent transport of c. Their trends are quite different 
in the two databases and provide a way to delineate between gradient and counter- 
gradient diffusion transport. As previously shown, the main differences between the 
two databases are the turbulence levels and the flame front wrinkling. 

6.2.1 Low turbulence level 

In the case of a low turbulence level and a low flame front wrinkling, (u.n) 5 
may be assumed almost constant and equal to its value at c = 0 where the normal 
component is equal to —1 in the mean propagation direction x: 

(u.n ) s = -Uf, (25) 

where Uf is the flow velocity at the flame front and may be estimated with the 
same argument as previously used by Bidaux & Bray (1994) to derive an expression 
for the turbulent flux of flame surface density (Eq. (13)): 

Uf = Z7 0 (1 + rc*)\ (26) 

where Uo is the fresh gases velocity and c * the flame reference level. Assuming a 
stationary flame brush and using mass conservation, an estimate of (u.n) 5 is: 

U 


{u.n) s = 


(1 + tc) 


(1 + TC*), 


(27) 



116 


A. Trouve, D. Veynante, K. N. C. Bray & T. Mantel 



Mean progress variable 


FIGURE 14. (u.n) s ( ), uncorrelated part (u) s .(n) s ( ), and correlated 

part (u.n) s — (u) 5 .{n) s ( ) plotted as a function of the mean progress variable, 

c, for the two available databases (top: CTR, bottom: Rutland). 

which is very well verified in the Rutland database. In this case, the flame front 
dynamics and the turbulent transport are mainly dominated by the mean flow 
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velocity and the thermal expansion. 

6.2.2 High turbulence level 

In the CTR database, the fluctuations at the flame front of the flow velocity and 
the flame normal vector may be assumed uncorrelated: 

(u”.n) s E = (u") s .(n) 5 S. (28) 

Combining (28) and (15), it may be easily shown that: 

„,(u".„) s E = (29) 

which suggests that (u".n) s might provide a way to analyze the turbulent diffusion 
transport of c. . 

The only unclosed term in Eq. (24) is (u") s . In the CTR DNS, the flame front 
movements are mainly dominated by the turbulent motions in the fresh gases (due 
to the higher viscosity, turbulent motions have a lower strength in the burnt gases). 
At the leading edge of the flame brush, the flame front is convected towards the 
fresh gases with a turbulent speed of the order of — u f where u f is the turbulent rms 
velocity in the fresh gases. At the trailing edge, the flame front is convected towards 
the burnt gases by the most energetic turbulent eddies inside the fresh gases with 
a turbulent speed of the order of u*. Then, a simple model for (u") 5 is: 

( u ")s = u ' (2c — 1) , (30) 

which is in agreement with the CTR database as shown in Fig. 15. 

From (28) and (30), one obtains: 

p b ( u".n) s E = -u'(2c- 1 )Vpc. (31) 

Using (29), assuming that c* = 0.5, the turbulent flux of c becomes: 

u n c" = —2uc(l — c), (32) 

which clearly corresponds to gradient turbulent transport. 

6.2.3 Mechanisms for turbulent transport diffusion 

From the previous analysis, the turbulent diffusion of c comes from two different 
mechanisms. The first one is the dynamics of the flame front itself, mainly due 
to thermal expansion. This phenomenon (clearly apparent in the term pcV.U of 
Eq. (24)) induces counter-gradient diffusions (CGD) as shown by the classical Bray- 
Moss-Libby expression (6). Due to thermal expansion and without any other effect, 
{/ft is greater than U u leading to CGD of c. 

In contrast, with a sufficiently high turbulence level, the flame front movement is 
dominated by the turbulent motions in the fresh gases as already seen in the CTR 
database. At the leading edge, the flame front is convected towards fresh gases, 
inducing negative velocities near the front, positive values of (u.n} 5 , and U b lower 
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FIGURE 15. Comparison of (u") s and the modeled expression (30) in the CTR 
database. : CTR database; : estimation from (30). 


than U u . At the trailing edge, intense turbulent motions in the^ fresh gases convect 
the flame front towards burnt gases and U u is greater than U t (see Fig. 16). In 
that case, the turbulent diffusion of c becomes gradient. 

To summarize, with a high turbulence level, the flame front dynamics is domi- 
nated by the turbulent motions and the mean progress variable is similar to any 
passive scalar, leading to a classical gradient diffusion turbulent transport. With 
a low turbulence level, the intrinsic dynamics of the flame is important, and the 
mean progress variable cannot be reduced to a passive scalar. A counter- gradient 
diffusion turbulent transport is then encountered. It is well known that counter- 
gradient diffusion is promoted by a high heat release (i.e. a high value of r). On the 
other hand, high turbulence levels induce gradient turbulent transport. The next 
step is to quantify the limits between gradient (GD) and counter-gradient (CGD) 
turbulent diffusion. 

6.3 A criterion for gradient/ counter- gradient turbulent diffusion 

The objective is to derive a criterion for turbulent diffusion transport of the 
Favre- averaged mean progress variable, c. Our analysis is based on Eq. (24). For 
GD, -V.pil"?' has to be positive at the leading edge of the flame and negative at 
the trailing edge. As previously shown, this term may be expressed as the sum of 
three contributions. At the trailing edge, the first term, ^(u^.n^E, is negative, 
the second one, {pb(w) s ~ ( P w )s J ) E vanishes and has a higher order, and the third 

one, pcV.U is positive. Then, a necessary condition for GD is: 
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FIGURE 16. Mechanisms of counter-gradient (left) and gradient (right) turbulent 
diffusion transport for c and S. 
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Ng = 


<^yL y_ s r(l +r) 
g6g r |(u".n) s | 


< 1 , 


( 38 ) 


evaluated for c = 1. This expression may be simplified by evaluating the turbu- 
lent flame speed, Sr, from the integral of the reaction rate over the flame brush. 
Assuming that (pw)s = one obtains: 


S T =^f-Pu(w) s . (39) 

Then the proposed criterion becomes: 

Ng = < i. (40) 

|(u".n) s | 

A similar analysis, conducted with Eq. (23), leads to a similar criterion: 


N'b = Ms 


1 + r(2 + r) 

l(u-n>sl 


< 1. 


(41) 


Table 1. Numerical estimates of the proposed criterions from the two available databases 


Database 

u 1 

SL 

T 

St 

(c= 1.) 

|(u" 

(c = 

•“)sl 

1-) 

Nb 

K 

CTR 

3. 

l. 

3. 

1.7 

8. 


3. 

4. 

2. 

Rutland 

1. 

l. 

2.3 

1.8 

3. 


1. 

7.6 

3.6 


Table 1 provides the values of the numbers Ng and Ng, for the two available 
databases, using the assumption ( w)s = S L- Both Ng and N' B are lower in the 
CTR database (exhibiting GD) them in the Rutland database (with CGD), a result 
that is consistent with our expectations. Nevertheless, the limiting value is not 
found to be equal to 1. This is mainly due to some assumptions made in the 
derivation. In particular, expressions (36) for the flame surface density £ and (37) 
for the gradient of the mean progress variable c are not well verified at the trailing 
edge of the flame front and only provide an order of magnitude. The criterion is 
also somewhat approximate because it corresponds only to a necessary condition to 
have a gradient diffusion transport. Two quantities are not exactly known in this 
criterion: the surface average of the flame front displacement speed, (w) s , which 
differs in general from the laminar flame speed s L (see Trouve & Poinsot 1994); and 
(u".n) 5 which probably involves wrinkling scales of the flame front. Nevertheless, 
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a rough estimate may be made assuming {w)s = &l and |(u".n) 5 | = w # . Then, the 
proposed criterion becomes: 




U' 

37 


(42) 


which corresponds to a horizontal line in the classical Borghi/Barrere diagram (fig- 
ure 3). Below this line, corresponding to low turbulence levels, turbulent transport 
occurs opposite to the direction of mean gradients. Above that line, turbulent trans- 
port occurs in the direction of mean gradients. This first simple criterion is also in 
agreement with the classical analysis of turbulent transport where counter-gradient 
diffusion is promoted by the increase of the heat release parameter r (Masuya 1986). 

The present analysis should be understood as a first attempt to differentiate be- 
tween gradient and counter-gradient turbulent diffusion phenomena in premixed 
flames. The criterion (42) is probably too simple, and the transition between 
counter-gradient and gradient transport as turbulence levels axe increased will not 
be abrupt and will involve a transition zone. It remains, however, that such a 
transition will probably occur as supported by our comparison of the two available 
DNS databases. Future work will be aimed at validating the present analysis and 
providing a more precise description of the transition between the two situations. 


7. Conclusions 

We use in this study direct numerical simulations to describe the coupling be- 
tween the transport equations for mean reaction progress variable, c, and flame 
surface density, E. We are particularly interested in the turbulent transport terms 
appearing in the equations for c and E since the modeling of those terms remains 
somewhat controversial: standard gradient transport approximations as proposed 
in the Coherent Flame Model, or additional transport equations for the turbulent 
fluxes as proposed in the Bray-Moss-Libby theory. Two different direct numeri- 
cal simulation databases, the CTR database and the Rutland database, have been 
used. Both databases correspond to statistically one- dimensional premixed flames 
in isotopic turbulent flow. The run parameters are significantly different and the 
two databases correspond to different turbulent combustion regimes (see Fig. 3). 

A systematic comparison between the different simulated flames reveals strik- 
ing differences in their turbulent transport properties. The Rutland case exhibits 
counter-gradient diffusion of c, a result that is consistent with the Bray-Moss-Libby 
theory. In contrast, the CTR database features turbulent c-fluxes that axe consis- 
tent with a gradient transport approximation. In addition, the turbulent diffusion 
of E is found to be strongly correlated to the turbulent diffusion of c. In the Rutland 
database, counter-gradient diffusion of c is occurring with counter-gradient diffusion 
of E, whereas in the CTR database, gradient diffusion of c and E is observed. A 
simple expression proposed by Bray to relate the turbulent flux of E to the tur- 
bulent flux of c is found to be valid (Eq. (15)). One important implication of this 
result is that the modeling of the turbulent flux of c cannot be made independently 
of that for the turbulent flux of E. 
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A detailed analysisof all the terms appearing in the conservation equation for the 
turbulent c-flux, ~pu n c n , is also performed. This analysis shows that while pressure 
effects act to promote counter- gradient diffusion, a result that is consistent with 
the Bray-Moss-Libby theory, they fail to prevail in the CTR case. Finally, a simple 
theoretical analysis based on a thin flame model is developed to explain the different 
turbulent transport properties observed in the CTR and the Rutland simulations. 
This simple theory distinguishes between situations where the flow field near the 
flame surface is dominated by the turbulent motions and situations where it is 
mainly determined by the dilatation occurring within the reaction zone. The CTR 
simulations belong to the first category, the Rutland simulation to the second. 
The first category of flames is expected to feature gradient turbulent transport, 
whereas counter-gradient phenomena will dominate in the second. The analysis 
also provides a criterion to determine the domain of occurrence of gradient /counter- 
gradient turbulent transport of c and S (Eq. (42)). This criterion suggests that 
counter-gradient diffusion is promoted by increasing values of the heat release factor, 
r, a result that is consistent with the Bray-Moss-Libby theory as well as by low 
turbulence levels, whereas gradient diffusion might prevail as turbulence intensities 
are increased. Future work will be aimed at validating and refining this criterion. 
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